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Abstract. The linear canonical transformations of geometric optics on two-dimensional 
screens form the group Sp(4, 5ft), whose maximal compact subgroup is the Fourier group 
U(2) F ; this includes isotropic and anisotropic Fourier transforms, screen rotations and gy- 
rations in the phase space of ray positions and optical momenta. Deforming classical optics 
into a Hamiltonian system whose positions and momenta range over a finite set of values, 
leads us to the finite oscillator model, which is ruled by the Lie algebra so(4). Two dis- 
Qh tinct subalgebra chains are used to model arrays of TV 2 points placed along Cartesian or 

polar (radius and angle) coordinates, thus realizing one case of separation in two discrete 
coordinates. The 7V 2 -vectors in this space are digital (pixellated) images on either of these 

S two grids, related by a unitary transformation. Here we examine the unitary action of the 
analogue Fourier group on such images, whose rotations are particularly visible. 
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O 
O 

\C> 1 Introduction 

O 

The real symplectic group Sp(4, 3?) of linear canonical transformations is widely used in paraxial 
• • geometric optics [1, Part 3] on two-dimensional screens, and in classical mechanics with quadratic 

.£h potentials. Rs maximal compact subgroup is the Fourier group U(2) p = U(1) F ® SU(2) p [2, 3]. 

These are classical Hamiltonian systems with two coordinates of position (q x , q y ) G 3ft 2 , and two 
coordinates of momentum (p x ,p y ) £ 3ft 2 , which form the phase space 3ft 4 , and which are subject to 
linear canonical transformations that preserve the volume element. The central subgroup U(1) F 
contains the fractional isotropic Fourier transforms, which rotate the planes (q x ,p x ) and (q y ,p y ) 
by a common angle. The complement group SU(2) F contains anisotropic Fourier transforms that 
rotate the planes (q x ,Px) and (q y ,p y ) by opposite angles; it also contains joint rotations of the 
(QxiQy) and (Px,Py) planes; and thirdly, gyrations, which 'cross-rotate' the planes (q x ,Py) and 
(q y ,p x ). The rest of the Sp(4, 3ft) transformations shear or squeeze phase space, as free flights 
and lenses, or harmonic oscillator potential jolts. 

This classical system can be quantized a la Schrddinger into paraxial wave or quantum 
models. Indeed, the group of canonical integral transforms was investigated early by Moshinsky 
and Quesne in quantum mechanics [4, 5], [6, Part 4] and by Collins in optics [7], and yielded 
the integral transform kernel that unitarily represents the (two- fold cover of the) group Sp(4, 3ft) 
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on the Hilbert space £ 2 (K 2 ). Yet it is the discrete version of this system which is of interest 
for its technological applications to sensing wavefields with CCD arrays. One line of research 
addresses the discretization of the integral kernel to a matrix and its computation using the fast 
FFT algorithm [8, 9, 10, 11]. This has the downside that the kernel matrices are in general not 
unitary, and thus no longer represent the group - non-compact groups can have only infinite- 
dimensional unitary representations [12]. 

We have developed a ' finite quantization 7 process to pass between classical quadratic Hamil- 
tonian systems, to systems whose position and momentum coordinates are discrete and finite 
[13, 14, 15, 16, 17]. The set of values a discrete and finite coordinate can have is the spectrum of 
the generator of a compact subalgebra of u(2), which is a deformation of the oscillator Lie algebra 
osc. When position space is two-dimensional, we use so(4) = su(2) x © su (2) y [17, 18, 19]. The 
purpose of this article is to elucidate the action of the Fourier group U(2) F on the iV 2 -dimensional 
representation spaces of the algebra so(4) that we use to realize pixellated images. 

Previously we have analyzed and synthesized pixellated images on sensor arrays that follow 
Cartesian or polar coordinates [18, 19, 20]; here we investigate the action of the Fourier group 
on finite pixellated images arranged along these two coordinate systems, understood as one 
example of separation of coordinates in discrete, finite spaces. In order to present a reason- 
ably self-contained review of our approach to two-dimensional finite systems, we must repeat 
developments that have appeared previously, which will be recognized by the citation of the 
relevant references. As research experience indicates however, each restatement of previous 
results yields a streamlined and better structured text, where the notation is unified and di- 
rected toward the economic statement of the solution to the problem at hand. In the present 
case, it is the action of the four-parameter Fourier group on Cartesian and polar-pixellated 
images. 

The classical Fourier subgroup of paraxial optics on two-dimensional screens [3] is described 
in Section 2. To finitely quantize this, we review the one-dimensional case, where u(2) is used 
to model the finite harmonic oscillator [13, 16, 17] in Section 3. The ascent to two-dimensional 
finite systems, Cartesian and polar, occupies the longer Section 4. There, under the aegis of so(4) 
we describe its finite position space separated in Cartesian [15] and in polar [18] coordinates, 
together with the two-parameter subgroups of the Fourier group that are within so(4) for each 
coordinate system. The new developments begin in Section 5, where we import the continuous 
group of rotations on images pixellated along Cartesian coordinates from their natural action on 
the polar pixellation. This serves to complete the action of the U(2) F Fourier group on Cartesian 
screens in Section 6. To translate the action of this group onto polar screens, we recall in Section 7 
the unitary transformation between Cartesian and polar-pixellated screens [19], thus finding the 
representation of the Fourier group on the circular grid as well. Finally, in Section 8 we offer 
some comments on the wider context in which we hope to place the separation of discrete 
coordinates in two dimensions. 

2 The classical Fourier group 

The classical oscillator is characterized by the Lie algebra osc of Poisson brackets between the 
oscillator Hamiltonian function h(q,p) := \{p 2 +q 2 ) and the phase space coordinates of position q 
and momentum p, 



{h,q} P = -p, 



{h,p} P = q, 



{q,p} P := 1, 




where 1 Poisson-commutes with all. The first two brackets are the geometric and dynamic 
Hamilton equations, while the last one actually defines the algebra of the system to be osc. 
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Figure 1. Action of the Fourier group on the sphere (from [21]). 



In two dimensions (q x ,q y ) £ and with the corresponding momenta (p x ,p y ) £ 3^ 2 , one can 
build ten independent quadratic functions 

QiQi', QiPi>, PiPi', i,i'e{x,y}, (2.2) 

that will close into the real symplectic Lie algebra sp(4, Four linear combinations of them 
generate transformations that include the fractional Fourier transform (FT) , which are of interest 
in optical image processing, 

isotropic FT F := \ (pi +p 2 y + q 2 x + q 2 ) =: \{h x + hy), 
anisotropic FT Fx := \ (p 2 x - p 2 y + q 2 x - q 2 ) =■ \{K - h y ), 
gyration F 2 := \{p x p y + q x q y ), 
rotation F 3 := \{q x p y - q y p x ). 

Their Poisson brackets close within the set, 

{Fx,F 2 } P = F 3 , and cyclically, 
{F ,F t } P = 0, i = 1,2,3, 

and characterize the Fourier Lie algebra u(2) F = u(l) F © su(2) F , with Fq central. This is the 
maximal compact subalgebra of (2.2), u(2) p C sp(4, K) [3]. In Fig. 1 we show the transformations 
generated by (2.4)-(2.6) on the linear space of the su(2) algebra, which leaves invariant the 
spheres F? + Ff + F| = F 2 . 

The Lie algebra u(2) and group U(2) will appear in several guises, so it is important to 
distinguish their 'physical' meaning. The classical Fourier group U(2) F of linear optics will be 
matched by a 'Fourier-Kravchuk' group U(2) K acting on the position space of finite 'sensor' 
arrays in the finite oscillator model. In one dimension, this model is based on the algebra u(2), 
whose generators are position, momentum, and energy; while in two, the algebra is su(2) x © 
su(2) = so(4). In all, the well-known properties of su(2) = so(3) will of use [22]. 



(2.3) 
(2.4) 
(2.5) 
(2.6) 



3 Finite quantization in one dimension 

The usual Schrddinger quantization (q^-Q = q-,p^-P = —\rjd/dq) of the phase space coor- 
dinates leads to the paraxial model of wave optics (for the reduced wavelength r\ = \/2ir) or to 
oscillator quantum mechanics (for rj = K). In one dimension, (2.2) yields the three generators of 
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the (double cover of the) group Sp(2, 3ft) of canonical integral transforms [4] acting on functions - 
continuous infinite signals - in the Hilbert space £ 2 (K). In two dimensions, the Schrodinger 
quantization of (2.2) yields the generators of Sp(4, 3?), and its Fourier subgroup is generated by 
the quadratic Schrodinger operators Fo, F\, F%, F% corresponding to (2.3)-(2.6), that we indi- 
cate with bars, and are up-to-second order differential operators acting on the Hilbert space. 
Finite quantization on the other hand, asigns self-adjoint N x N matrices q H > Q, p H > P to the 
coordinates of phase space. 

Since the oscillator algebra is noncompact, it cannot have a faithful representation by finite 
self-adjoint matrices [12]; we must thus deform the four-generator osc into a compact algebra, of 
which there is a single choice: u(2). We should keep the geometry and dynamics of the harmonic 
oscillator contained in the two Hamilton equations in (2.1), with commutators i[- , •] in place 
of the Poisson brackets {•, -} P , and in place of the oscillator Hamiltonian h, a matrix K added 
with some constant times the unit matrix. We call K the pseudo-Hamiltonian. The third - and 
fundamental - commutator [Q, P] will then give back K. In [13] we proposed the assignments 
of self- adjoint N x N matrices given by the well-known irreducible representations of the Lie 
algebra su(2) of spin j, thus of dimension 2V = 2j + 1, where the matrix representing position 
is chosen diagonal. Using the generic notation { Ji, J%, J3} for generators of su(2), the matrices 
are 

position: q \-t Q = J\, 

Qm,m> = mS m>ml , m,m' G {-j, -j + (3.1) 
momentum: p 1— > P = J2, 

P m ,m' = -ilVti -m)(j + m + l)5 m+ i,m' + ' l \ V U + m)(j-m + l)<5 m -i, m / , (3.2) 
pseudo-energy: K = J3, energy: h\-t H := K + (j + 

K m ,m' = \ Vti ~ m )U + m + l)<5 m+ i, m ' + \ \[Q + m)(j - m + l)<5 m _i, m '. (3.3) 
The commutation relations of these matrices are 

[K, Q] = -iP, [K, P] = iQ, [Q, P] = -iK, 

i.e., [Ji, J2] = — i/3, and cyclically. (3.4) 

The central generator 1 of osc is corresponded with the unit N x N matrix that generates a U(l) 
central group of multiplication by overall phases. This completes the algebra u(2) realized by 
matrices acting on the linear space C N of complex iV-vectors. 

The spectra E of the matrices Q of position {q}, P of momentum {w}, and K of pseudo- 
energy {k} are thus discrete and finite, 

S(Q) = E(P) = Tl(K) = {-j, -j + 1, . . . , j}, 

with j integer or half-integer. We define the mode number as n := K + j € {0, 1, . . . , 2j}, and 
energy is n + \. The Lie group generated by {1;Q,P, K} is the U(2) group of linear unitary 
transformations of C N . Finally, the sum of squares of these matrices on C N is the Casimir 
invariant 

C:= Q 2 + P 2 + K 2 =j(j + 1)1, 

In Dirac notation, C N is a spin-j u(2) representation, where the eigenstates of position q and of 
mode n = k + j satisfy 

Q\q)i = q\q\, K \ n ) 3 = ( n ~ j)\ n ) 3 , q\-j, n\l 3 . (3.5) 
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The position eigenvectors \q\ form a Kronecker basis for C N where the components of |n) 3 will 
be provided, from (3.3), by a three-term relation, i.e., a difference equation that is satisfied by 
the Wigner little- d functions [22] for k = n — j. (The reader may be disconcerted for having J\ 
diagonal, whereas almost universally J3 is declared the diagonal matrix; formulas obtained by 
permuting 1 1— > 2 1— >■ 3 1— > 1 are unchanged.) The overlaps between the two bases in (3.5) are the 
finite oscillator wavefunctions, 

^n(q) := 1 (q\n) 3 = 4-^(1^) = *W n ~ j) ( 3 - 6 ) 
(-If i(2j\ ( 2j 




2j y V n ) \q + j 



K n (q + j;l,2j). (3.7) 



This overlap is expressed here [23] in terms of the square root of a binomial coefficient in q, 
which is a discrete version of the Gaussian, and a symmetric Kravchuk polynomial of degree n, 
K n {q) = K n (q; \,N — 1) = K q (n — j) [24]. We have called the $ n (g)'s Kravchuk functions [13]; 
they form a real, orthonormal and complete basis for C N . The Lie exponential of the self-adjoint 
matrix K generates the unitary U(l) group of fractional Fourier-Kravchuk transforms, whose 
realization we shall indicate by U(1) K . 

For future use we give the general expression of the Wigner little-d functions for angles 
/3 € [0, 7r] as a trigonometric polynomial [12, 22] 

C'.rJ/ 3 ) = VU + m')\{j - m')\(j + m)!(j - m)\ 

^ (-l) m '- m+fc (cos i/3)^'+ m - m '- 2fc (sin ^p)m'-m+2k ^ 



= 4n>JrP) = (-l) m - m 'd^ mj - m/ (p), (3.9) 

where due to the denominator factorials, the summation extends over the integer range of 
max(0, m — m!) < k < min(j — m' , j + m); for m — m! > the reflection formulas (3.9) apply. 



4 Two-dimensional systems 

The two-dimensional classical oscillator algebra 0SC2 is sometimes taken to consist of the Poisson 
brackets between 1, q x , q y , p x , p y , and h x , h y ; and sometimes an angular momentum m = 
q\Vi ~ Q2P1 is included, Poisson-commuting with the total Hamiltonian h = h x + h y . The finite 
quantization of two-dimensional systems deforms 0SC2 to the Lie algebra su(2) © su(2) = so(4) 
in both cases. (The subalgebra of 1 need not concern us here.) To work comfortably with the 
six generators of so(4), we consider the customary realization of J^i = —Ji\i by self-adjoint 
operators that generate rotations in the (i, i') planes, which obey the commutation relations 

[Ji,i', Jk,k'] = ^i',kJi,k' + 0~i,k>Ji>,k + 0~k,iJk>,i' + 0~k' ,i> Jk,i) , 

that we summarize with the following pattern: 



Jl,2 


Jl,3 


Jl,i 






^2,4 






J3,A 



(4.1) 



A generator J^y has non-zero commutator with all those in its row i and its column i' (reflected 
across the i = i' line); all its other commutators are zero. The asignments of position, momentum 
and pseudo-energy with the su(2) = so(3) matrices (3.1)-(3.3) is 
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K 


-P 




Jl,2 


Jl,3 




h 


-J2 




Q 






^2,3 






Jl 



4.1 The Cartesian coordinate system 

Passing from one to two dimensions can be achieved prima facie by building the direct sum 
algebra su(2) x . © su(2)^, which is accidentally equal to so(4). This isomorphism is shown in 
patterns by 



V 















Qx 1 






Qx H 




P'x 















(4.2) 



where the square super-index of all generators, X?, indicates the identification between the 
so(4) generators Ja> with the Cartesian coordinates and observables. 

Since the ^-generators commute with the y-generators in (4.2), using (3.5) a Cartesian basis 
of positions and also a basis of modes can be simply defined as direct products \q x ,qy)\ '■= 
l&LkvXy and \ n x,n y )^ := \n x ).Jn y ) 3y , 



Qx \Qxi Qy)i 



Qy \lx, Qy)i 



K x \nx,n y )^ = (n x - j) \n x ,n y )^, 



K° \n x ,n 



y/3 



Qy \Qx: q y )i 
= (% - 3) \n x ,n y )s, 



(4.3) 
(4.4) 



where the pseudo-energy eigenvalues are n% — j = K %\-j- The two-dimensional finite oscillator 
wavefunctions in Cartesian coordinates are thus given as a product of two x I / n (g , )'s from (3.6), 
(3.7), 



*^,n u fe,%) = a i(qx,q y \n x ,n y }^ = tfn«(g*)*nv(«v)» 



(4.5) 



on positions q x , q y \_j and of mode numbers n x , The positions can be accomodated in the 

square pattern of Fig. 2 a, and the modes in the rhombus pattern of Fig. 26. The Cartesian 
eigenstates of the finite oscillator are shown in Fig. 3. 

Since K x and K y generate independent rotations in the (Q X ,P X ) an d in the (Q y ,P y ) planes 



respectively, their sum K := K x + K y = J\p, transforms phase space isotropically. Thus we 
identify K G so(4) as the generator of a group U(1) K of isotropic fractional Fourier-Kravchuk 



C n2 , and corresponding with the operator 2Fq G U(2) f in (2.3). That 



transforms in C N x C N 

sum commutes with A := K x — K y = J34 G so(4), which generates skew-symmetric Fourier 
rotations by opposite angles in the x- and y- phase planes; so A corresponds with 2F\ G u(2) F 
in (2.4). However, a counterpart of the 2F3 generator of rotations in the (q x ,q y ) and (p x ,p y ) 
planes cannot be found within so(4), and neither can the gyrations in (2.5). These will be 
imported in the next section. 



4.2 The polar coordinate system 

The six generators of the Lie algebra so(4) can be identified with positions and modes following 
an asignment different from the Cartesian direct sum of the previous section. We indicate the 
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Figure 2. Cartesian eigenvalues following the asignments of position and energy generators in (4.2) 
within the algebra so(4). (a): Position eigenvalues (q x ,q y ) in (4.3). (b): sum and difference of mode 
eigenvalues (n x ,n y ) in (4.4). 





Figure 3. Basis of x- and y-mode eigenstates on the Cartesian array, VPS^n,, (q x , %) in (4.5). 
new generators with a circle super- index, X°; as in the classical case, a generator of rotations M 



between the x- and y-axes should satisfy the commutation relations 

[M,Q°]=iQ°, [M,Q° y ] = -iQ° x , (4.6) 
[M, P°] = iP°, [M, P°] = -iP°, (4.7) 

while the isotropic Fourier generator K = should rotate between position and momentum 
operators, 

{K,Q° x ]=iP°, [K,P°]=-iQ x , (4.8) 

[K,Q° y ]=iP°, [K,P°] = -iQ° y , and (4.9) 

[K,M}=0. (4.10) 

The commutator (4.10) asserts that K and M can be used to define a basis for C n2 with the 



quantum numbers of mode and angular momentum, that will be given below. Expressed in the 
pattern (4.1), a new so(4) generator assignment that fulfills these requirements is 
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— P° 


-p° 




Ql 


Q°y 






M 



(4.11) 



This assignment satisfies the conditions (4.6)— (4.10), but implies the further commutators 

[Q x ,P°]=iK=[Q° y ,P°], (4.12) 
[Q x ,Q y ]=iM=[P°,P°], (4.13) 
[Ql,P°]=0=[Q o y ,P°]. 



Of these, (4.12) echoes the su(2) nonstandard commutator in (3.4), while the commutator (4.13) 
is also nonstandard, and indicates that Q° and Q° cannot be simultaneously diagonalized. 

First, we find operators with quantum numbers corresponding to radius and angle; we use 
the subalgebra chain 



K 


-K 






Ql 


Ql 






M 



Ql 


Ql 




M 



D 



M 



(4.14) 



When both su(2)'s in (4.2) have the same Casimir eigenvalue j(j + 1), the principal Casimir 
invariant is Y ^tv = 2j (j ' + 1)1. In so(4) there is a second invariant, Y Ja'Jkk'i which 

i<V ' (i<i')^{k<k>) 

is identically zero in these 'square' cases [22]. Let us now consider the Casimir operator of the 
subalgebra so(3) C so(4) in (4.14), which is 



R(R + l) := (Q° x y + (Q a y y + M 



(4.15) 



The Gel'fand-Tsetlin branching rules [25] determine that C n2 then decomposes into subspaces C p 
that are irreducible under this so(3), where (4.15) exhibits the eigenvalues [22] 



p€{o,i,...,i-i}. 



Although R is not an element of the algebra so(4), we shall identify it as the radius operator. 
The final link in the reduction (4.14) is M, whose eigenvalues in the ^-representation of so(3) are 

m G {— p, —p + 1, . . . , p}. The total number of distinct eigenvalue pairs (p, m) is Yl + 1) = 

(2j + l) 2 = iV 2 , the same as for the N x N square grid, and shown in Fig. 4a. Thus we define 
the radius- angular momentum (ra) eigenvectors 



R\p>™)tU. = p\Pi*$U> M\p,m)° A = m\p,m)° A , p\ 3 , 



ml 



(4.16) 



and we adopt p as the radial position coordinate. 

In the last step we use the discrete Fourier transform matrix, F p in each (2p+ l)-dimensional 
subspace, to pass between angular momenta and angles, and thus build the Kronecker basis of 
states 4>) & localized at a definite radius p and the 2/3 + 1 equidistant angles 4>k as 



\P, 9k) 



— y 

j2p + l ^ 

v ' m=—p 



exp(-im0 fc ) \p,m)° A , 



(4.17) 
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(Jb ' 


p o > 
• i 


p 

• 


• < 






• • 

• 
• 


m • ^ — 
• 

• 



Figure 4. Polar eigenvalues following the asignments of radius and angle in (4.17) within the algebra 
reduction so(4) D so(3) D so(2) in (4.14). Left: Eigenvalues (p, m) in (4.16). Right: Positions of radius p 
and angles 4>k according to (4.17) aligned by ijj p = 0. 



1 P 

Ip,™)ra = /2 1 Yl ex P(+ im ^fc) \P, 
k=—p 

for fc := 2vr/c/(2p + 1) + -p<k< p, 

where the ipp's are fixed but arbitrary phases. In Fig. 4b we show the resulting arrangement 
of iV 2 points (p,4>k) thus defined. Note that the Fourier transformation in (4.17) is linear and 
unitary but is not an element of the group SO(4): it has been imported to act on each of the 
2p + 1-dimensional so(3) irreducible representation subspaces p G {0, 1, ... , 2j}. 

Having the position Kronecker eigenstate basis \p, (j>k) e f° r polar coordinates, we now define 
the basis of mode and angular momentum MA, eigenbasis of the commuting operators K = 
Kx + Ky with total mode number (pseudo-energy k = n — 2j), and of M = Kx + Ky with 
angular momentum m, placed in the pattern (4.11). We build these eigenstates \n, m)^ A asking 
for 

We observe that whereas the RA states in (4.16) are classified by the Gel'fand-Tsetlin chain of 
subalgebras so(4) D so(3) D so(2), the MA states in (4.18) follow the chain so(4) = su(2)©su(2) D 
u(l)©u(l), which in ordinary quantum theory entails the coupling of two spin-j representations 
to total spin p. The overlaps between the MA and RA states should therefore be Clebsch- 
Gordan coefficients ~ C mi, m 2 , m, with mi = ^(k + m) and m2 = —\{k — m), adding to the 
total m [22]. In the present construction though, the subalgebra chain (4.14) reduces along the 
'lower' subalgebras, and this differs from the original Gel'fand-Tsetlin reduction that reduces 
along the upper ones; also, generally one counts su(2) multiplet states from the top down, and 
we have ordered them from the bottom up. The overlap of RA and MA states entails an extra 
phase that must be computed carefully. It is [18, 19] 

R ° A (p, m\n + j, m)° MA = <p(j, p, k, m) C [ m+K)/ ^ ( m _ K)A P m , 
tp(j, p, k, m) := (— iy +p exp[i^7r(K + \m\ — m)\. 

The overlap between the states of mode n = k + j and angular momentum m with the polar 
Kronecker basis of radius p and angle (j)k yields the discrete polar oscillator wavefunctions 



^°n,m{P^k) ■= Q {pAk\n- j,m)° MA 

2== £ e-^j, p, k, m) C { m+n _ j)/2t J m _ n+i)/2> p m . (4.19) 



v ' m=—p 
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o 

## 

I- 




n 



m 



Figure 5. Basis of mode-angular momentum cigenstates on the polar array of radius and angle, 
K, m (P,<M in (4.19). 



These states are accommodated in a rhombus (n,m), similar but distinct from the rhombus in 
Fig. 2, which was classified by (n x ,n y ); here it consists of the eigenvalue pairs 



ii 



ii 



2j 
' 



m G {— n, —n + 2, . . . , n}, 

m £ {—4:j + n, — 4j + n + 2, . . . , 4j — n}. 



(4.20) 



In Fig. 5 we show the C n2 basis of MA states. Note that due to the Clebsch-Gordan selection 
rules, states of a given angular momentum m are nonzero only at radii p > \m\. 

The generator K in (4.11) generates rotations between both position and momentum opera- 
tors, and corresponds to twice the isotropic Fourier transform generator 2Fq in (2.3). The action 
of K.{oS) := exp(— koK) on the MA basis of C N is thus 



fC(u) |n,m)°, 



e -2i(n-2,v | n ,m)° A . 



(4.21) 



Similarly, rotations are generated by angular momentum, 7Z(9) := exp(-WM); and since M 
is twice the generator F3 E SU(2) F in (2.6), the vectors in the polar basis (4.19) of C N are 
multiplied by a phase with the double angle, 



K{6) |n,m)° y 



e- 2ime |n,m)° A . 



Under these rotations, images f(p, 
transform into 



(4.22) 

Jp,(f>k\f) on the polar-pixellated screen will thus 



fe(p, ■= f(p, fa + 0)= Q (p, 4> k \n{6) 1 /) 

= Jfi, fa\ exp(-i^M) |p, m)° A ° A (p, m\f) 
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Figure 6. Rotation of images on the polar screen. 



R°(p; 4>k, fa;0)f(p, 4> k ') = f{p, 4>k-0), 

k'=-p 



where for each radius pig' 7 there is a [2p + 1) x (2p + 1) matrix representing the same rotation 



>ki 9k' 



Jp, fap, m)° A exp(-i6>m) ° A (p, m\p, fa) 



2p 



1 P 

— ^ exp[-im(6 - 4> k + fa)] 

m=—p 

1 sin[(p+ \){9-<p k + fa)} 



2p+l sini( 



^k + fa 



(4.23) 



These are circulating matrices, functions of 4>k — 4>k' = 27r (A; — k')/{2p + 1) modulo 2tt, and 
periodic in k, k' modulo 2p + 1. For each radius, the 'unit' rotation angle is 9 = 2ir/{2p + 1), 
and for multiples I thereof, the matrix (4.23) is nonzero at the diagonal k = k' + I. In Fig. 6 we 
give an example of such rotation. 

Isotropic Fourier transformations and rotations on the polar screen are produced by genera- 
tors within the so(4) algebra in the pattern (4.11). However, the pattern also informs us that 
with linear combinations of K and M, we can not gyrate the planes (Q°,P°) and (Qy, P£) jointly 
as with 2F2 in (2.5); also missing is the anisotropic Fourier transform generated by 2Fi in (2.4), 
which was natural in the Cartesian basis. These transformations will be imported on C N in 
Section 7. 



5 Importation of rotations on the Cartesian screen 

We noted above that in the subalgebra chain (4.14), the generators of isotropic Fourier trans- 
formations and rotations, K 2Fo and M 2F3 are domestic to so(4), while anisotropic 
Fourier transformations and gyrations, corresponding to 2F2, are foreign. Now, in the Cartesian 
subalgebra decomposition of so(4) in (4.2), the two independent Fourier transform generators 
(2.3), (2.4), Kx Fq + Fi and Ky f> Fq — Fi are domestic to so(4), while those of gyrations and 
rotations, 2F2 and 2F3, are foreign. Since we cannot complete a fully domestic U(2) K Fourier- 
Kravchuk group in correspondence with the Fourier group U(2) F C Sp(4, K), we must import 
the missing transformations. Such importation was used in (4.17) with the (2p + 1) x {2p + 1) 
discrete Fourier transform matrix. 

We now build the group of rotations on the Cartesian grid by importing SU(2) transformations 
[22, Chapter 3] from the continuous model. Rotations of an image should respect the energy of 
each formant mode n = n x + n y = k + 2j in the \n x , %) 3 D basis of C N , and transforms real 
images into real ones, while mixing states with different eigenvalues p := \{n x — n y ), Ai|"(5 2 j 
horizontally across the rhombus in Fig. 2b. The proposed imported action of ^3 in (2.6) on the 
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Cartesian modes stems from (3.3), replacing j \— > ^{rix + n y ) = and m = fi \— > -^{jix — n y ), 
namely 

M\n x ,n y )° = \Jn y (n x + 1) \n x + l,n y - 1) 3 D + \jn x {n y + 1) \n x -l,n y + 1) 3 D . (5.1) 

We must pay attention to the fact that F3 in (2.6) is one-half of the angular momentum operator 
QxPy — QyPx that is the generator of finite rotations, TZ(0) exp(— 219F3). On the Cartesian 
mode basis states its action involves the standard Wigner little-d function (3.8) through 

K(0) K, n y )° = d Tx-n y )/2,(n' x -n> y )/2( 2d ) K' n 'yV • 

K+ n 'y= n 

This is a real linear combination of the Cartesian mode basis where n and \i are bound within 

the rhombus of Fig. 2b, for integer ri | q- 7 , H™{f/ 2 * n the l° wer half and for n\^-, ^|^y 2 -j 2 m ^ ne 
upper one. 

In Fig. 7 we show rotation of a Cartesian-pixellated image f(q x ,q y ) = i{lx,Qy\f) to 

fe(<h»q v ) ■= a i(q x ,q y \nO) I/) = q y \n{6) \n x ,n y )° %n x ,n y \q' x ,q' y )°%q' x ,q' y \f) 
= R D (Qx,qy;q' x ,q' y ,0)f(q' x ,q l y ). 

The subgroup of rotations 7£(#) G SU(2) K is thus represented by the iV 2 x N 2 matrices 
R a (Qx , q y ; q' x , q' y ; 0) : = a -lq x , q y \ n x , n y ) 3 %i x , % |72.(0) I n' x , n^) 3 D , | q£ , q y ) ° 

= E^ >% fe,%)<;(2^)*°, )< (^,^) ! (5.2) 

where /i = ^(n x — n y ) and ^' = ^(n^. — n' y ) are bound by n x + n y = n = n' x + n' y , and belong to 
the same row in the rhombus of Fig. 2b. 



6 Completion of U(2) K on the Cartesian screen 

Having imported a unitary representation of the group of rotations onto pixellated images on 
the Cartesian screen, and having the domestic group of fractional Fourier transforms, we can 
complete the Fourier-Kravchuk group U(2) K on this screen. From (2.4) and (4.2), F 1 O \{K° - 
K y ), the Lie exponential A(4>) -B- exp(— 2i(pFi) acts on the Cartesian kets (4.4) through phases, 

A(4>)\n x ,n y )° = exp[-2i(j)(n x — n y )]\n x ,riy)i . (6.1) 

Hence, for images f(q x ,q y ), 

U(qx,q y ) ■= A{(j)) f{q x ,q y ) = ^ A a (q x ,q y ;q' x ,q' y ;(j))f(q' x ,q' y ), 
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with the matrix kernel 

A a (q x ,q y ;q' x ,q' y ;<p) = ^n x ,n y fe> %) exp[-2i(f){n x - %)]*^,n a {.<&, ( 6 - 2 ) 

where the sum over /i = \{n x — n y ) preserves n. 

Now, having two of the three generators of SU(2) K , we can produce the third: gyrations 
G(ip) O exp(— 211PF2), through 

On images f(q x ,Qy), Fourier-Kravchuk gyrations will act through a matrix kernel, 
U{qx,q y ) ■■= G(ip) f(q x ,q y ) = ^ G a (q x ,q y ;q x ,q y ;i(;)f(q x ,q y ), 

CP (q x ,q y ; q' x , q' y , = £ ^ fe, q y )e~^ d n J* (2V0e +i ^' /4 *°, , n , ^ ^ 

where, as in (5.2) and (6.2), the sums over /i = \(n x — n y ) and fi' = \(n' x — n' y ) preserve n. In 
continuum optics, gyrations acting on the Hermite-Gauss beams transform them into Laguerre- 
Gauss ones of the same mode number [21, Fig. 4]. 

Finally, the isotropic K,(oj) <-> exp(— 2iuFq) in (4.21) is domestic to the so(4) algebra (4.2), 
and completes the U(2) K group with 

fUqx,q y ) ■= f(q x , q y ) = Yl ^ D fe>^;^^;w)/(<^.,<^), 

q'x, q'y 

K D (q x ,q y ;q x ,q y ;uj) = ^ ^n x ,n y (qx, q y ) exp[-2iu(n x + n y )}^° xjny (q' x , q' y ). 

The elements of the Fourier group U(2) F = U(1) F ® SU(2) P , where the factors are comple- 
mentary, are customarily parametrized by Euler angles as 

V(uj;4>,9,ip) := exp(— iuFo) exp(— i^Fs) exp(— 16F2) exp(— itpFs), (6.3) 

and its matrix elements between eigenstates (i,fj.\ and |fc, //), with eigenvalues /i, fi' under F3 
and 1 under Fq, the latter being the irreducible representation label are the well-known Wigner 
Big-D functions 

As we indicated in Section 3, by permuting 1 1— > 2 1- > 3 1- ?■ 1 in (6.3) and using (4.21), (4.22), 
and (6.1), we can write the elements of the isomorphic Fourier-Kravchuk group U(2) K as a pro- 
duct of Fourier-Kravchuk transforms and rotations, 

v(u;<f>,e,ii)) = K{\uj)A{\(j))n{\e)A{\ii)). 

Its matrix elements between the Cartesian mode eigenstates \n x ,n y )P of K x and KyVfiW. be 
thus, 

>,, % |P(u;;0,0^) \n' x ,n' y }° = e'^' 2 ^ D^_ nM< _ n , )/2 (0, 9, j>), (6.4) 
with the total mode n = n x + n y = n' x + n' as before. 
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The domestic and imported transformations properly mesh, and we see that the Fourier- 
Kravchuk group U(2) K is indeed represented unitarily and faithfully by (6.4) on C N . Its action on 
images f(q x , q y ) pixellated on the Cartesian screen can be found from here, writing VL := (0, 6, ip), 
through a real similarity transformation by the matrix formed with the C N basis (4.5) of 
Cartesian mode Kravchuk functions, 

fu,Q.{<lx,qy) ■= T>(u},n) f(q x ,q y ) = ^ D a (q x ,q y ;q' x ,q' y ;uj,n) f{q' x ,q' y ), (6.5) 

where, with n = n x + n y = n' x + n' y and u\q the kernel is 
D a (q x ,q y ;q x ,q y ;uJ,ty 

Tlx yfly 

n' x n' y 

and these N 2 x N 2 matrices are unitary representations of U(2) on the Cartesian grid of points 
in Fig. 2b. With this result we now turn to the polar screen. 

7 The Fourier group on polar screens 

In the subalgebra chain of so(4) that produces the polar screen (4.11), the commuting rotations 
and isotropic Fourier-Kravchuk transforms are domestic. To complete the Fourier-Kravchuk 
group U(2) K on the polar screen, we must import either gyrations or the anisotropic transform. 
Since these transformations have been realized already in the Cartesian screen basis, we should 
find a unitary map between both screens. This transformation has been studied in [18, 19], 
and consists in identifying the eigenstates of mode and angular momentum in both bases; in the 
polar basis these are \p, m)^ A in (4.16) shown in Fig. 5. In the Cartesian basis such states will be 
constructed now as eigenvectors of 1Z{9) = exp(— \9M) with eigenvalues e~ l9m , or equivalently 
of with eigenvalues |m, corresponding to the eigenvalues of |M in (2.6). These Cartesian 
states will be linear combinations, respecting total mode number n = n x + n y , of all states 

\ n i m )MA := Cn'™n y \n x , %)i D , 

n x +n y =n 

M \n, m)° A = m \n, m)° A , K \n, m)° A = \n |n, m)° A . 
From (5.1), the coefficients obey the difference equation 

y/ny(n x + l)C^i >ns _i - \mCn x m n y + y/n x (ny + 1)^^1,^+1 = °> 

which is the ubiquitous su(2) three-term recursion relation of the Wigner little-d functions for 
angle \i\ (around the 1-axis [22]), now with j O \n = ^(n x + n y ) and \m ■H- \{n x — n y ). We 
thus define MA states having angular momentum m on the square grid as 

|n,m)° A := £ ^^^^p^) K, n y )?. (7.1) 

n x +n y =n 

These states can be accomodated in a rhombus (n,m), exactly as in (4.20). 

From (7.1) follows the definition of the Cartesian basis of mode and angular momentum MA 
eigenstates 

^n, m {<lx,qy) ■= a i(qx,q y \n,m)^ A (7.2) 
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Figure 8. Basis of mode-angular momentum eigenstates on the Cartesian array, An, m (Qx, %) m (7-1). 

E>K-%)/4//2 (±Tr\^ n (n a \ 

e fl m/2,K-»j)/2^ 1 ^ ni '"»W 1 ' ? W 

n x +n y =n 

= A*_ m (q x , qy y = (-l)n+m A n m{ _ qx _ qy) = (-l)^+^A°_„ i _ m fe )% ). 

This transforms the Cartesian Kravchuk functions into what we called Laguerre-Kravchuk func- 
tions in [21, Fig. 4]. The basis of states (7.2) is shown in Fig. 8, which can now be identified 
with the basis - also of mode and angular momentum - in Fig. 5. We can now import the 
equivalence between the MA bases 

\n,m)^ A = \n,m)° lA . (7.3) 

A pixellated image on the Cartesian screen, f a (q x ,q y ) £ C n2 , can be thus unitarily trans- 
formed into an image on the polar screen, f (p,<j>k) G C n2 through the transformation 

fo(p,4>k) = e {p,4>k\f) 

= XI Q(P^k\qx,q y )? a i(qx,q y \f) = ^2 U (P^k;qx,q y )f a (qx,q y ), (7.4) 
where the transform kernel, using (7.3), is 

U(p,(j)k;q x ,q y ) ■■= Q {p,<f>k\qx,q y }° 

= Y^ &(p,<Pk\n,m)° MA °iA(n,m\q x ,q y )° =J2^°n, m (pAk)^, m (qx,q y T- 

n,m n,m 

The transformation of images inverse to (7.4), from the polar to the Cartesian screen, is 

fa(.Qx,q y ) = ^V{q x ,q y ;p,4> k )f {p,(j) k ), 

p,<i>k 
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Figure 9. Unitary map (7.4) of an image (the letter R) on the Cartesian screen onto the polar screen. 
V(q x ,q y ;p,(j) k ) = ^ A° m {q x , q y )^ m (p, </> k )* = U(p,<p k ;q x ,q y )* . 

n,m 

We note that these transformations map real functions onto real functions. In Fig. 9 we show 
an example of the map (7.4) on the 0-and-l image of the letter R. 

The action of the Fourier-Kravchuk group U(2) K on images on the polar screen can be now 
derived from its action (6.5), (6.6) on the Cartesian screen: 

fw,si{p,4>k) ■= fi) f(p, <j) k ) = De (Pi4>k;p,4>k';u,n)f(p,4> k <), 

where the kernel is a representation of U(2) K on the polar screen 

D Q (p,c/) k ;p',(j) k r,u;,n) = ^ U(p,(j) k ;q x ,q y )D a (q x ,q y ;q' x ,q' y ;uj,n)V(p ,(j) k >;q x ,q y ). 

qx,q y 

Q'x,<ly 

We have thus far not found a more compact analytic form for the representation of the Fourier 
group on the discrete coordinates of radius and angle. 

The iV 2 x N 2 matrix U = (V) ' that is needed above, is arduous to calculate but its numerical 
values need be computed only once for each size N of the screen, and will serve to transform any 
image from the Cartesian to the polar screen and back. For N = 32 there are ~ 10 6 values to 
be stored, and the transformation of an image as that in Fig. 9 between the screens involves this 
number of sums and products. Commercial image-rotating algorithms store the original image 
and interpolate the few Cartesian pixels around each geometrically determined point. Certainly, 
the Fourier-Kravchuk group of transformations of images presented here does not provide a fast 
algorithm, but it is unitary, and hence does not loose information, as interpolation invariably 
does; it is exact. 

8 Concluding remarks 

Willard Miller Jr. is recognized as having initiated the modern study of systems whose governing 
equation allows for separation in various coordinate systems using Lie algebraic reduction meth- 
ods [26]. These equations are differential because space has been considered continuous. Finite 
spaces with discrete coordinates, used as foundation for finite Hamiltonian systems character- 
ized by governing difference equations, have not been considered - to our knowledge - in the 
same context. The present paper reviews our work in the particular case of a two-dimensional 
Hamiltonian system whose two governing Hamilton equations identify as a finite harmonic os- 
cillator, mothered by the compact algebra so(4). And since it happens that algebra has two 
distinct subalgebra reductions, we can relate them with Cartesian an polar coordinates. 

The difference equations whose solutions are the finite oscillator wavefunctions and the 
Clebsch-Gordan coefficients are contained in so(3) and so(4), and entail special Kravchuk and 
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Hahn polynomials. Other finite polynomials, such as those of Meixner and Pollaczek, also 
appear in related harmonic oscillator models [27], and thus could be interpreted in terms of 
two-dimensional physical and optical models. Further, finite three-dimensional coordinate sys- 
tems beyond Cartesian and circular cylindric ones may be of interest in applications. But we 
should end with a word of caution regarding applications: the two-dimensional arrays that 
purportedly measure the quality of Laguerre-Gauss (polar) laser beams generally use sampled 
Laguerre-Gauss functions; the question of whether these, or the mode-angular momentum finite 
functions, are the most efficient to find the mode and angular momenta of the actual beams, 
seems to favor the sampled functions [28, 29], although they do not form orthonormal bases for 
the space of pixellated images. But only the discrete function bases synthesize all images exactly. 

The continuous two-dimensional harmonic oscillator is a superintegrable system, which is 
known to separate generally in elliptic coordinates, of which the Cartesian and polar are limits. 
We have been unable so far to find a corresponding finite array that would allow for pixellation 
following ellipses and hyperbolas, although there are three enticing leads: diagonalization of 
the so(3) operator J 2 + aJ^ [26]; a line of functions with a parameter which joins continuously 
Kronecker deltas with Clebsch-Gordan coefficients [30]; and the production by holographic 
means of laser beams whose nodes follow elliptic coordinates [31]. There are also 'geometrical' 
reasons to doubt that this is possible, however: whereas in Cartesian coordinates all pixel sizes 
are equal, and in polar ones they are almost so (except very near to the center), trying to 
visualize these two pixellations as limits of a finite elliptic pixellation presents some problems 
that we have not yet been able to overcome [32]. 
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